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, It is evident that models of the knee should match the observational phenomenology. In this talk I 

discuss a few aspects of phenomenology, which are important not only for the understanding of the 
knee origin, but also for the general problem of the origin of cosmic rays. Among them are the shape of 
the energy spectrum, its irregularity, the sharpness of the knee and its fine structure. The classification 
, of models is given and some examples of the most recent models are discussed. The most probable 

conclusion deduced from this examination is that the knee has an astrophysical origin and the so called 
'source' models of the knee are most likely among them. 
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The total number of models pretending to explain the origin of the knee exceeds twenty 
and the number of people who propose these models is about thirty. The number of ideas 
r% , is definitely larger, but I deliberately distinguish ideas and models, because sometimes the 

ideas are not developed at all to be considered as a model. The quality of the development 
is also quite different for different models. To the great extent it is due to the lack of astro- 
physical or particle physics data which could constrain the models, but there are also some 
evident subjective reasons. Certainly it is our duty to examine as much consequences of the 
proposed idea as we can and to some extent it is our fault that we ignore some evident ex- 
perimental or observational facts which contradict the proposed idea. There is a tradition in 
our P.N. Lebedev Physical Institute that if you propose the idea you should also propose the 
test which can kill your idea. We must admit that there is not a lot of authors who propose 
such tests. 

In this talk I examine some features which should be explained by models of the knee 
and requirements which should be satisfied. 

2. Phenomenology 

Usually we consider just three main characteristics of cosmic rays (CR) which should be 
predicted by the model and checked by the comparison with experimental data - they are 
the energy spectrum, the mass composition and the anisotropy. Among them there might 
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be some phenomenological sub-features, which are sometimes crucially important. Quite 
a good survey of models with their regard to the observed features of the energy spectrum 
and the mass composition is given in. 1 Here I focus only on the problems related to the CR 
energy spectrum. 

Usually we examine such features of the spectrum as the absolute intensity and slopes 
below and above the knee. I should like to draw the attention to some other sub-features of 
the spectrum shape viz. its negligible irregularity, then the sharpness and the fine structure 
of the knee. 

2.1. The irregularity 

Since the discovery of extensive air showers (EAS) it has been found that their size spec- 
trum has the power law character in wide energy ranges below and above the knee. The 
theory of shock acceleration explained this power law and since then the educated guess 
was that it is supernova explosions which create the shocks, which in turn accelerate par- 
ticles and produce a regular power law of the CR energy spectrum. However, we should 
remark that a stochastic character of supernova explosions in space and time prevents the 
formation of the regular spectrum. The subsequent propagation of CR in the interstellar 
medium (ISM) only emphasizes the irregularity since it results in the dominant contribu- 
tion of nearby and recent supernova remnants (SNR) to CR observed at Earth. The accurate 
Monte-Carlo realization of the scenario, where only Galactic SNR give birth to observed 
CR, shows the extremely irregular character of the spectrum with the very low probability 
to observe the good power law (Figure la). 

We consider the irregularity as a very important phenomenological feature and propose 
as the measure of the irregularity the parameter T = <j\ og j, which is the standard deviation 
from the mean intensity taken in a logarithmic scale. As an example, Figure lb shows 
the irregularity as the function of primary energy for the scenario of standard supernova 
exploding sporadically in our Galaxy and producing different spectra, shown in Figure 
la. It is seen that due to the decreasing lifetime of higher energy CR in the Galaxy the 
smaller number of recent SNR contributes to the intensity and the irregularity of the spectra 
increases. 

2.2. The sharpness 

The remarkable sharpness of the knee has been mentioned even in the first publication 
by Kulikov and Khristiansen in 1958. 2 They have written: The presented data indicate 
with a great probability the possible sharp change of the spectrum in the studied range 
of shower sizes .... In 3 we suggested to use as the measure of sharpness the parameter 
S = — gn^jTATp ' where I is the shower intensity and N is the EAS size. The advantage of 
this second derivative is that it depends neither on the absolute intensity nor on the slopes of 
the spectrum below and above the knee. The world survey of 40 spectra 4 showed that inspite 
of the big spread of sharpness values nearly all of them exceed the value 0.3 expected 
for the original Galactic Modulation Model (Figure 2). This important phenomenological 
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Fig. 1 . (a) A sample of 50 CR energy spectra produced by stochastic explosions of supernovae and hypernovae 
in our Galaxy. Each spectrum is formed by explosions of 49000 standard supernovae and 1000 hypernovae. (b) 
The mean spectrum and its irregularity (standard deviation of the CR intensity). The rise of the irregularity with 
the energy is clearly seen. 

characteristic should be taken into account in every reasonable model of the knee and never 
ignored. 



2.3. The fine structure 

Fluctuations in the EAS development smear out the irregularity of the primary CR energy 
spectrum, but not down the regular power law with T = 0. Due to the stochastic nature 
of the CR sources there must be 'structure'af some level. The well known examples of 
such structures are the widely discussed knee and the ankle. The search for other structures 
is difficult since the deviations from the power law are evdently small. One must have a 
good statistics to find the structure in the single experiment. To use the data accumulated in 
many experiments one has to eliminate their difference in the background and calibration. 
We satisfied the first requirement using the deviation from the running mean. To satisfy the 
second requirement we referred the position of the deviation to the position of the knee 
in each individual experiment. The position of the knee has been determined as the point 
of the maximum sharpness in the spectrum. This approach is absolutely legitimate if you 
study the shape of the spectrum. 

We processed 40 EAS size spectra and 5 Cherenkov light spectra available at that time 
and the result is shown in Figure 3. Besides the apparent excess at the knee there is a second 
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dotted line: SSM, cr exp = 0.0 
dashed line: SSM, a = 0.1 
full line: SSM, C7 exp = P 0.2 
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Fig. 2. The sharpness of the knee in the Cherenkov light (the dashed box) and EAS size spectrum vs. the 
atmospheric depth of the measurements. The origin of symbols is given in 4 . The lines are the predictions of 
the Single Source Model with an account for the errors of the measurements u logNs : - dotted line, 0.1 - 
dashed line, 0.2 - full line. The horizontal dashed line at S = 0.3 indicates the prediction of the classic Galactic 
Modulation Model (GMM). It is seen that the bulk of measured sharpness values exceeds the predictions of the 
GMM. 

excess ('peak') at \og(E/E knee ) « log(iV e /iV e fc " ee ) « 0.6 seen both in the Cherenkov 
light and in the EAS electron size spectrum. We consider this result as the strong evidence 
for the fine structure of the knee, which should be taken into account by any reasonable 
model of the knee. 

The origin of this peak will be discussed later. However it is clear that such structure is 
nearly impossible to get in the classic Galactic Modulation Model with its smooth leakage 
of CR from the Galaxy. 



3. Classification of models 

All models of the knee can be tentatively divided into three groups: astrophysical models, 
interaction models and phenomenological models. The astrophysical models in turn can be 
subdivided into source models and propagation models. 

The models of the first group explain the origin of the knee by astrophysical reasons: 
either by the change or the evolution of CR sources, or by the change of the propagation 
and the leakage from the Galaxy. The interaction models claim that there is no sharp change 
of the primary CR energy spectrum. They refer the steepening of the observed EAS size 
spectrum and the change of their characteristics to the change of the particle interactions at 
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Fig. 3. The excess over the running mean: (a) in the primary energy spectrum expected in the SSM; (b) in the 
Cherenkov light spectra; (c) in the EAS size spectra. The lines indicate the predictions of the SSM for the versions 
with O, Fe (full line) and He, O (dashed line) responsible for the first (knee) and the second peak respectively. 



PeV energies. The phenomenological models don't specify new sources or new processes 
above the knee, but try to explain the steepening by the existing features of the primary 
CR or an EAS development. Sometimes the association of the model with any group is 
ambiguous, and this classification is given just for the orientation. Below I'll give some 
examples of the models, which were discussed in recent years. Certainly the list of these 
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examples is not complete and I ask those authors who will not be mentioned to apologize 
me for my ignorance. 

3.1. Astrophysical models 

3.1.1. Source models 

The most developed conception is that sources of CR below the knee are supernova (SN) 
explosions. The shock waves produced by the explosions can accelerate protons up to max- 
imum energies E max w 0. 1 - 0.4 Pe V. The knee is due to the impossibility for SN to ac- 
celerate particles up to higher energies. This cut-off at E max is rather sharp and by this 
way one can explain the visible sharpness of the knee. However the immediate questions 
to these models are: 

(i) are SN standard enough to give about the same E max to give the sharp knee ? 

(ii) what is the origin of CR above the knee? 

Different source models give different answers to these questions. I should like to speak 
about our Single Source Model later and here just mention the model of P.L.Biermann. 5 He 
allows that there might be just a few ordinary SN remnants (SNR) which give dominant 
contribution at the knee and ensure the formation of the sharp knee. He even assumes that 
their explosions due to special conditions created during the collapse can be standard. 6 
The cosmic rays above the knee are produced by another type of SN, which progenitors 
are more massive, hot and emit powerful particle wind before the explosion. The typical 
representative of such type of sources are well known Wolf-Raye stars. The powerful shock 
wave propagating through the bubble created by the hot wind accelerates particle up to EeV 
energies. 

Similar idea is developed by L.G.Sveshnikova, 7 who examined the variety of SN types 
and concluded that the dominant contribution to CR below the knee is given by just SNI 6c 
type and above the knee - by explosions of more massive and hot stars which she called 
'hypernovae' though this name is often associated with SN explosions accompanied by 
gamma-ray bursts (GRB). 

The interesting model of the cosmic ray acceleration by blobs of relativistic plasma is 
proposed by A.Dar, R.Plaga and A.De Rujula. 8-13 Its main subject is the origin of gamma 
ray bursts, but there are important points in this scenario which have relevance to our dis- 
cussion. The sources of cosmic rays in this model are 'cannonballs' - relativistic plasmoids 
which are emitted by the core collapsed supernovae and stopped in the Galactic Halo or in 
the Intergalactic Medium (IGM). By this way authors avoid difficulties connected with the 
need to introduce different sources below and above the knee - the sources are the same 
and the knee is connected with the increasing number of reflections from the cannonball, 
needed to accelerate particles to higher energies. Since the model of the knee is not suf- 
ficiently developed it is difficult to say whether it explains the sharpness of the knee or 
not. As for the irregularity the problem seems not to be severe since CR sources can be 
even more numerous than SNR since each of them can emit a few cannonballs. There is no 
problem with the isotropy and the absence or presence of the GZK cutoff, since the sources 
are mainly distributed in the Galactic Halo 9 or in the IGM, 12 but not in the Galactic Disk. 
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There are other models which claim that there is no need to introduce another kind of 
sources beyond the knee. SNR are able to accelerate particles up to hundreds PeV due to 
compression and magnification of magnetic fields inside the shell. 14 However these strong 
fields decrease with time as soon as the SNR expands. SNR turned out to accelerate parti- 
cles to higher energies during the shorter time at the very beginning of their history. It gives 
rise to the steepening of the energy spectrum and imitating the knee. 15 ' 16 

There are so called compact source models which are the combination of the source 
and interaction models. They examine the possibility for protons emitted by the source to 
lose the energy in photopion reactions and for nuclei - to disintegrate after the interaction 
with intense radiation fields around these compact sources. 17 ' 18 Though the knee formed by 
these process is sharper than in the propagation models considered below (S w 0.35 instead 
of 0.3) it is not sufficient to explain the observations. Moreover these models predict the 
extremely light mass composition beyond the knee, which contradicts the observations. 

I think all the mentioned source models have problems when compared with the shower 
phenomenology. Models which require different sources below and above the knee usually 
propose more energetic and less abundant sources of higher energy CR above the knee. 
They inevitably meet problems with the irregularity. The small number of high energy 
sources or shorter time of the acceleration to higher energies increase the sensitivity of 
the cosmic ray energy spectrum to the particular space-time distribution of the sources 
and reduce the possibility to form the regular power-law spectrum. Source models need 
the additional mechanism (Halo or IGM) which could smooth the irregularity of the CR 
spectrum below and above the knee. On the other hand the models with SNR accelerating 
particles up to sub-EeV and EeV energies give just the smooth steepening of the CR energy 
spectrum and are not able to explain the sharpness of the knee. 

3.1.2. Propagation models 

Propagation models were the first suggested after the discovery of the knee. 19 They as- 
sociate the steeper spectrum above the knee with the increased leakage of CR particles 
from the Galaxy when their gyroradius becomes comparable with the transverse size of the 
Galactic Disk. Latest developments of this idea included models with introduction of the 
Hall diffusion 20 22 and an anomalous diffusion in the ISM of quasi-fractal type. 23 

All they have problems with the sharpness of the observed spectrum. Original mod- 
els had the sharpness of the constituent spectra not greater than 0.6. For the mixed mass 
composition with the abundance of constituent nuclei equal to their abundance at TeV en- 
ergies the model sharpness decreased down to 0.3. In more sophisticated models with the 
Hall diffusion and a special geometry of regular magnetic field in the Galactic Disk and 
Halo the authors 20 21 created the small intensity bump in the knee region and by this way 
increased the sharpness of constituent spectra up to 1.0. The composite spectrum with a 
mixed composition will also have a greater sharpness but no more than 0.7. It is much less 
than the required value of about 3. 

The second problem is the observed fine structure of the spectrum. If the propagation 
models predict the smooth constituent spectra with a continuous steepening, it is impossible 
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to construct a total spectrum with discontinuities which are seen as the fine structure. 

3.2. Interaction models 

The interaction models introduce new processes in the interaction of CR particles at PeV 
energies to justify the sharp steepening of the EAS size spectrum. Some of these models re- 
quire the new physics either to create new heavy particles with large cross-sections, which 
decay into yet unobservable high energy muons and neutrinos, 24 ' 25 some required emission 
of technihadrons, gravitons 26 or very many low energy particles. 27 The model with an en- 
hanced production of charmed particles 28 slows down the development of the atmospheric 
cascade so that the bulk of its energy is dissipated not in the atmosphere but reaches the 
mountain or even the sea level and is absorbed in the earth. 

The experimental confirmation of neutrino oscillations and therefore an existence of 
the neutrino mass inspired the creation of the models where cosmic ray particles interact 
with relic neutrinos in the gravitationally attracted neutrino halos surrounding galaxies and 
galaxy clusters. 29 By assuming the neutrino mass in the region of eV it is easy to obtain the 
energy threshold for the inelastic interaction in the PeV region. 

In my mind the drawback of these models is not only that most of them require the new 
physics in the energy region very close to that already studied in Fermilab, but that they 
often don't specify the new process and don't develop their interaction model to incorpo- 
rate and test it with well developed simulation codes such as CORSIKA, like it has been 
done for more conservative models (QGSJET, SYBILL etc.). These models can explain the 
existence of the knee, sometimes even its sharpness by introducing a rapid rise of the cross- 
section above the threshold, but never other phenomenology of EAS accumulated hitherto. 
In this respect I should better attribute nearly all of them to ideas and not to models. 

3.3. Phenomenological models 

Phenomenological models of the knee just extrapolate the observed phenomenological 
characteristics of CR or EAS into the high energy region and demonstrate that with cer- 
tain approximations they give the steepening of the energy spectrum which they associate 
with the knee. 

To such type of models I attribute the model of Tsallis et al. 30 They applied the formulae 
of generalized non-extensive statistics to the observed CR energy spectrum and have found 
that it is possible to fit its behaviour with the bend at sub-GeV and the steepening at PeV 
energies in the total range of 34 decades of the differential intensity. The knee in this model 
appears as a crossover between two fractal thermal regimes. 

To the phenomenological class of models belongs also the multi-component model of 
Biermann and Ter-Antonyan. 31 They obtain parameters of the constituent nuclei spectra 
by multivariate fitting of the electron and muon size spectra measured by KASCADE and 
MAKET ANI experiments. 

Another phenomenological model is proposed by Stenkin, 32 who connects the steep- 
ening of the EAS size spectrum with the non-linearity of the logiV e (log.Eo) dependence 
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which in turn is connected with the change of the EAS structure as soon as the shower max- 
imum approaches the observation level. If N e oc E s where s is the decreasing function of 
Eq then even for the power law primary energy spectrum I(Eq) oc E^ 1 with a constant 7 
the slope index a = 2 of the EAS size spectrum I{N e ) oc N~ a increases with N e and the 
spectrum steepens imitating the knee even for the constant slope of the primary spectrum. 

I would attribute to phenomenological models also the so-called poly-gonato model of 
J.R.Horandel. 33 He noticed that spectra of constituent nuclei in primary CR become flatter 
with an increasing primary mass. He extrapolates these spectra up to their rigidity cutoff 
and found that it is possible to describe not only the knee but the general behaviour of the 
primary energy spectrum up to sub-EeV energies. At these high energies the substantial 
fraction of cosmic rays can be trans-iron elements. 

Since these models are purely phenomenological and not pretend to explain the origin 
of cosmic rays they have no problems with the spectrum irregularities. However they may 
have problems with the observed fine structure of the spectrum and definitely with the 
sharpness of the spectrum if they will not assume the sharp cutoff of the constituent spectra. 

4. Status of the Single Source Model 

It would be unfair if I'll tell nothing about the present status of the Single Source Model 
(SSM). Me and Prof. A.W.Wolfendale proposed it in 1997 3 and updated in 2001 . 4 The main 
idea was that the knee is due to the contribution of just one single, recent and nearby SNR. 
Since it is nearby and recent, CR from this SNR have a flat energy spectrum with a sharp 
cutoff. 34 At PeV energies their intensity becomes comparable with the total CR background 
from all other distant and old SNR and their contribution creates the bump which forms 
the knee. Due to the sharp cutoff of the spectrum beyong its maximum rigidity R max 
the knee is sharp and different constituent nuclei which have different maximum energy 
Emax = ZR max give rise to the fine structure of the spectrum. 

Actually the energy spectrum of CR created by the Single Source is the result of the fit 
of theoretical spectra 34 to the existing experimental data from 40 independent EAS electron 
size spectra, 9 muon size spectra, 3 hadron size and energy spectra and 5 Cherenkov light 
spectra. Comparing this spectrum with the spectrum of CR expected from SNR of different 
ages and placed at different distances from the Earth we estimated the most likely intervals 
of distances and ages of the Single Source. 35 They appeared to be rather narrow around 300 
pc and 100 kyear. Soon after that american astronomers using the triangulation technique 
have found that the pulsar B0656+14 associated with the SNR Monogem Ring is situated 
exactly at the indicated range of distances and ages as required for our Single Source. 36 

In 35 we presented arguments why the nearby SNR is not seen in gamma rays. It is 
because being nearby it is situated in our Local Superbubble with its low density of ISM 
which is about 3 • 10~ 3 cm~ 3 . Moreover it is not a discrete, but an extended source with 
the angular radius of 20° and the problem of the correct subtraction of the background 
from such a big area is not trivial. Many experiments have already looked for an excessive 
intensity of CR from this region and found nothing. 37 " 46 Nevertheless the situation is con- 
troversial, since recently the Moscow University 47 and Tien-Shan 48 groups revealed weak 
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Fig. 4. The structure of the EAS size spectra beyond the knee (full line). Dashed line shows the EAS muon size 
spectrum measured in the Moscow University 50 

EAS excesses in the Monogem Ring area. If there are indeed CR excesses at sub-PeV and 
PeV energies from the Monogem Ring region they can evidence for the activity of the pul- 
sar B0656+14 which emitted high energy CR soon after its birth and which were confined 
within the SNR for a long time of about 80 kyear. 49 High energy gamma-quanta which 
create these excesses are produced by CR on the local fluctuations of the interstellar gas 
density. 

The energy spectum of CR produced by the pulsar B0656+14 has the sharp peak at the 
rigidity of 0.25 PV which is very close to the R ma x ~ 0.4 PV, required by the SN model of 
Berezhko et al 34 and can contribute to the observed knee sharpness, but its contribution to 
the CR intensity at the knee does not exceed 15%. However, this pulsar can be responsible 
for CR beyond the knee at energies above 10 PeV. Recently we analysed the region beyond 
the knee and have found an extended excess of EAS intensity in the broad region slightly 
above 10 8 GeV 50 (Figure 4). The existence of this excess gives an argument in favor of 
the assumption that the pulsar B0656+14 contributes to the CR intensity above 10 7 GeV. 
So that there is an impression that the SNR Monogem Ring with its associated pulsar 
B0656+14 can be a serious contender for the Single Source responsible for the PeV CR 
intensity bump and the knee itself. 

It is not easy to kill the SSM, since its parameters were taken to fit the existing phe- 
nomenology. It is possible, however, to check and to kill the initial mass composition pre- 
dicted by SSM. Originally we attributed intensity peaks at 3 and 12 PeV to O and Fe. If 
we are right then we cannot expect another peaks beyond 12 PeV. However if the observed 
peaks are actually due to He and O, then we can expect the existence of the sharp Fe 
peak at 40 PeV. So far there is no indication of this peak in the experimental data from 
KASCADE. We look forward at the new data from KASCADE-Grande. 
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5. Summary 

In the summary I stress again that the models which pretend to explain the origin of the knee 
should explain not only the high isotropy, the tendency of the primary mass composition to 
become heavier above the knee, but also such phenomenological features as the sharpness, 
the fine structure of the knee and the relatively small irregularity of the energy spectrum 
both far below and far above the knee. At the moment none of the models can satisfy these 
requirements. 

Source models which require less abundant and more energetic sources have problems 
with the large irregularity of the spectrum in particular above the knee. The same difficulty 
have the models with relatively abundant sources like pulsars which emit high energy CR 
for a short time after their creation. Definitely we have to incorporate the object or the 
mechanism which level out irregularities originated from the stochastic space-time distri- 
butions of sources. The first candidates for such an object could be the Galactic Halo or 
IGM. 

All propagation models have not the sufficient sharpness of the knee even for the con- 
stituent spectra. As the consequence they should have problems with the observed fine 
structure of the spectrum. They take for granted the regular power law spectra with no ex- 
planation of the origin of this regularity, therefore they can also have hidden problems with 
the irregularity of the spectrum. 

The phenomenological models don't pretend to explain the origin of CR and the knee 
- they use different input spectra and find the best fit to the observed phenomenology of 
the showers. The sharpness, the fine structure depend on the shape of these input spectra so 
that there is no actual explanation of the observed phenomenology. 

The interaction models, besides that they require 'the new physics', are not sufficiently 
developed to be tested by the phenomenology. It is clear that modifying the threshold be- 
haviour of the new process it is possible to fit both the sharpness and the fine structure of 
the spectrum, but the value of this exercise is not very high. Since the main subject of in- 
teraction models is the knee they do not explain the small irregularity of the CR spectrum 
outside the knee region. 

Schematically the possibility for different models to match the phenomenology of EAS 
in the knee region is outlined in Table I. This table shows that the astrophysical models 



Table 1 . Correspondence between models and the phenomenology 



Model/Phenomenology 


Sharpness 


Fine structure 


Irregularity 


Source 


+ 


+ 


? 


Propagation 






not examined 


Phenomenological 


+ 


+ 


not examined 


Interaction 


+ 


+ 


not examined 



are the most developed and some of them can match the observed phenomenology. Among 
astrophysical models source models are the most plausible. For nearly all models we need 
to introduce a mechanism which could smear out the irregularity of the energy spectrum 
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due to the stochastic nature of the sources, viz. Galactic Halo, Intergalactic Medium or so. 
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